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We consider a ferromagnetic spin-1 Bose-Einstein condensate (BEC) dispersively coupled to a 
unidirectional ring cavity. We show that the ability of a cavity to modify, in a highly nonlinear 
fashion, matter-wave phase shifts adds a new dimension to the study of spinor condensates both 
within and beyond the single-mode approximation. In addition to demonstrating strong matter- 
wave bistability as in our earlier publication [L. Zhou et at, Phys. Rev. Lett. 103, 160403 (2009)], 
we show that the interplay between atomic and cavity fields can greatly enrich both the physics of 
critical slowing down in spin mixing dynamics and the physics of spin-domain formation in spinor 
condensates. 
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I. INTRODUCTION 

Experimental realization of spinor Bose-Einstein con- 
densates (BEC) has opened up a new research direc- 
tion of cold atom physics [l| , in which superfluidity and 
magnetism are simultaneously realized. Compared to 
scalar condensates, spinor condensates possess unique 
features: (i) The spin-dependent collision interactions al- 
low for the population exchange among hyperfine spin 
states; (ii) The spinor condensate is described by an or- 
der parameter with vector character and therefore may 
exhibit spontaneous magnetic ordering. These give rise 
to spin-dependent phenomena such as coherent spin mix- 
ing, spin textures and vortices, spin waves and spin do- 
mains. These phenomena have been extensively studied 
in theory |9l421| and demonstrated by a few pioneering 
experimental works [H-Q. 

In the study of spinor BEC, it has been found that 
magnetic field plays an important role, particularly via 
the quadratic Zeeman effect. Coherent control of the 
spin-dependent behavior has been achieved by tuning 
magnetic field. These include the control of the oscil- 
lation period and amplitude of coherent spin mixing &- 
la. [ll| - KL4| . formation of spin domain structure [ll-la. ll7H2l| 
and quantum phase transitions between different mag- 
netically ordered states [I, H, [l3| • 

In another frontier of cold atom research, recent exper- 
imental progress have realized strong coup ling of BEC to 
electromagnetic modes of optical cavity [22, [23|- This 
heralds a new regime of cavity quantum electrodynam- 
ics, where a cavity field at the level of a single photon 
can significantly affect the collective motion of the atomic 
samples, hence opening up new possibilities in manipu- 
lating ultracold atomic gases with cavity-mediated non- 
linear interaction. Previous works focused on the in- 
terplay between the cavity field and the atomic exter- 
nal degrees of freedom — the center-of-mass motion of 
scalar condensates |24| - [33 |. The ground state and collec- 



tive excitations [27|, |28| , cavity induced Mott insulator- 
superfluid phase transition [3l| and cavity optomechanics 
[33 | were theoretically investigated in detail. Such a sys- 
tem was also shown to have the potential applications in 
probing atomic quantum statistics in optical lattices and 
atomic quantum state preparation [33|. Experimentally, 
optical bistability at few-photon level has been observed, 
which is made possible by the strong atom-photon cou- 
pling 0,111. 

In our recent work [SJI, a system of a spin-1 BEC 
trapped inside a unidirectional ring cavity was studied, 
where the cavity couples to the atomic internal spin de- 
grees of freedom. We examined the equilibrium proper- 
ties of this system under the single-mode approximation 
(SMA) and showed that the interplay between the atomic 
spin mixing and the cavity light field can lead to strong 
matter-wave and optical bistability simultaneously. Our 
current work is an extension of Ref. [34i] . Here we will 
conduct a more complete investigation by including the 
study on the non-equilibrium properties and the collec- 
tive excitations of the system. We will also examine the 
validity of the SMA and show that, when SMA becomes 
invalid, spatial domain structure will form in the spinor 
condensate. This study will help us gain insight into such 
properties as the spinor dynamics, dynamical stability, 
spin domain formation, etc. 

The rest of the paper is organized as follows. Section 
II introduces the theoretical model. Section HI is de- 
voted to a discussion of spinor dynamics under the SMA, 
where both equilibrium and non-equilibrium properties 
are studied. The validity of the SMA is examined in Sec. 
IV by investigating the modulational stability of a ho- 
mogeneous system. We then present results showing the 
formation of spin domain structure in the ground state 
in the regime where the SMA becomes invalid. Finally 
we conclude in Sec. V. 



II. MODEL 



The equations of motion then read 




FIG. 1: (color online) Schematic diagram showing the system un- 
der consideration. An 7^ = 1 spinor condensate is trapped inside 
the cavity using an optical dipole trap. The population of different 
spin components can exchange via spin mixing. The cavity is co- 
herently driven by an external laser with amplitude Ep and decays 
with a rate K. The cavity field is vr polarized and is dispcrsively 
coupled to the atomic system. 



ihip±^ £ + Uo\a\ + C2 {n± + no - n^) i^± + C2fpQtp^, 

(la) 

ihij^a = [C + C2 (n+ + n_)] -i/jo + 2c2'0+V'-V'o, (lb) 

a^[iSc-iUo{N++N^)- K]a + ep, (Ic) 



where C = p^ /2ma + Vt (r) 4 
part of the Hamihonian, n - 
atomic density, (5c = Wp — oJc 



cqu is the spin-independent 
: n+ + riQ + n_ is the total 
is the cavity detuning rela- 



tive to the external laser field, and Uq = g / {oJp — LOg) is 
the effective atom-photon coupling, with g being atom- 
cavity mode coupling constant. Further, since the cav- 
ity decay rate k is typically much larger than the spin 
oscillation frequency, in what follows, we adiabatically 
eliminate a from Eq. (jlc|) . replacing a in Eq. (Ilap with 



a{t) 



K-i[5c-Uti {N+ + iV_ 



(2) 



Figure [T] is a schematic of our model in which a spinor 
BEC with hyperfine spin Fg = 1 confined in an opti- 
cal dipole trap by a spin-independent trapping poten- 
tial Vt (r) is placed inside a unidirectional ring cavity. 
The cavity is driven by a coherent laser field with am- 
plitude Ep and frequency Wp, far detuned away from 
the Fg = 1 ■'r^ F^ — \ transition frequency uja such 
that the atom-photon interaction is essentially of dis- 
persive nature. Further, the cavity is assumed to sup- 
port a single 7r-polarized electromagnetic mode charac- 
terized by a frequency Wc, a decay rate k. The selection 



rule then allows states 



Fg = l,m, 



= ±1) to be coupled 



to the corresponding states in the excited manifold with 
the same magnetic quantum numbers \Fi, = 1, vrie = ±1), 
while forbids state \Fg — 1, nig — 0) to make dipole tran- 
sitions to any excited states as shown in Fig. [1] How- 
ever, population in state \Fg = l,mg = 0) can be dis- 
tributed to other sublevels via the two-body s-wave spin 
exchange collisions. The total s-wave interaction are de- 
scribed by two numbers, cq = Anh? (2a2 + ao) /Snia and 
C2 — Anti^ (a2 — ao) /^nia, where ttIq is the atom mass, 
and a/ is the s-wave scattering lengths in the hyperfine 
channel with a total spin / = or 2. Since the anti- 
ferromagnetic case of ^'^Na has been considered in our 
previous work [3J|, here we will focus on the ferromag- 
netic case of ®^Rb spin-1 BEC where C2 < 0. 

In this work, we take the standard mean-field ap- 
proach, describing the cavity field with a complex am- 
plitude a (t) which amounts to assuming the cavity field 
to be represented by a coherent state, and the spinor 
condensate with the order parameters ip- (r, t) , ipo (r, t) , 
and ^+ (r, t) , which represent the wavefunctions in mag- 
netic sublevels nig = —1,0, and -1-1, respectively. This 
treatment is justified when the condensate atom number 
Na = J Ua (r) dr (where Ua = {ipal is the atom num- 
ber density) in magnetic sublevel a are sufficiently large. 



One may immediately observe from Eqs. ([T|) that the 
dispersive interaction between cavity photons and the 
condensate atoms introduces an effective quadratic Zee- 
man energy shift, [/o|ap, to nig = ±1 states relative to 
the TUg = state. However, unlike the Zeeman shift due 
to an external magnetic field or to a strong off-resonant 
laser field 16], a key feature of this effective shift is that 
it is sensitive to the spin population distribution of the 
condensate, as manifested by Eq. ([2]). As such, it gener- 
ates a new effective spin-dependent interaction which in 
turn induces a new set of nonlinear phenomena in spinor 
condensate. In what follows, we will describe in detail 
such new phenomena. 



III. SPIN DYNAMICS UNDER SMA 

In this section, we consider the spin dynamics under 
the assumption of SMA. This describes, for example, a 
condensate whose size is smalle r than the spin healing 
length ^s defined as ^s = h/ y^2ma \c2\n which repre- 
sents a length scale over which a local perturbation in 
spin density gets forgotten. Under the SMA, each spin 
component shares the same spatial wavefuntion (j) (r) ac- 
cording to 



-0Q (r, t) = VN(j) (r) VP^exp [-i {^it + Oa)] , q = ±, 0, 

(3) 
where 9^ is the phase, pa is the population normalized 
with respect to the total atom number TV = Y]^ Ng, 
and (f> (r) is the solution to the time-independent Gross- 
Pitaevskii equation: C(f> — ficj), where fi is the chemical 
potential and (f) (r) satisfy the normalization condition 
/dr|0(r)|^ = l. 

By inserting Eq. ([3]) into Eqs. (|lap and (jlbp . we arrive 



at a set of equations 



^ = 2AaPo-\/(l ^ /°o)^ -m^sinf 



dpo 
dr 
dd 
d^ 



= -2- 



t/o|a| 



2AaX 



, (1 - po) (1 - 2po) - m^ _ ^ 
1 — 2/9o H , — cos t 



^{l~ Pof --m? 



(4a) 



(4b) 



which describe the dynamics of a mixed state in which 
none of the spin component vanishes, where 9 — 29q — 
9^ — 0_ is the relative phase, m ^ p+ — p- the magneti- 
zation, and T = Kt the dimensionless time. In Eqs. Q, 
we have also introduced other dimensionless quantities 
given by 



Aa = 



Nc2jdr\^{r)f NUo Sp . S^ 
, Uo ^ , ?7 = — , dc = — ■ 

K K K K 



To facilitate our study below, we follow Refs. [l2|, |35| 
and use dpo/dr = —2dH/d9 and d9/dT = 2dH/dpo to 
construct, in terms of two conjugate variables po and 9, 
the following mean-field Hamiltonian H 



H = XaPO 

where 



1 - Po + Y (1 - paf -m'^cos( 



Uipo), 
(5) 



U {po) = — arctan [C/q (1 - po) - 5^ 
represents the cavity-mediated atom-atom interaction. 



A. Equilibrium Property: Bistability 

In this subsection, we will use Eqs. ^ to study the 
equilibrium property of a condensate in the parameter 
regime that supports bistability. As can be seen from 
Eq. (|4ap . at steady state, there are two branches of sta- 
tionary solutions: one with 9 = Q (the in-phase state) 
and the other with 9 = tt (the out-of-phase state). The 
in-phase state always has a lower energy for C2 < and 
we will therefore only focus on the in-phase state in this 
work. In addition, we will restrict ourselves to the case 
with zero magnetization m = 0, i.e., we only consider 
the case where there are equal number of rUg — 1 and 
rUg = — 1 atoms. 

Under these conditions, the intracavity photon num- 
ber can be found, by combing the stationary solution of 
Eq. (|4b[) with Eq. ^, to obey the following transcenden- 
tal equation 
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FIG. 2: (color online) (a) Mean intracavity photon number JQJ 
and (b) the normalized spin-0 population po versus cavity-pump 
detuning 5c for a steady-state solution with 5 = 0. The ones 
represented by the red dashed lines correspond to dynamically un- 
stable solutions. The vertical dotted lines indicate the position of 
the first-order transition which occurs at Sc = —4.75. 



where A = Uq/2 - Sc and x = U^/4:NXa- It is well- 
known that when 77^ \x\ > 8v3/9, the system will display 
bistable behavior [36]. 

Figure [2ja) shows how the intracavity photon number 
changes with detuning Sc, based on a set of realistic pa- 
rameters: Xa = -6.8 X 10"^ [37], Uo = -5, 77^ = 5, and 
A^ = 2 X 10^. With this set of parameters, rj'^ \x\ is found 
to be around 2.3, which is above the threshold value 
8^3/9 « 1.54. Indeed, for -4.9 < Sc < -4.6, the sys- 
tem supports three stationary solutions. The dynamical 
properties of these solutions can be studied with the stan- 
dard linear stability analysis. Substituting po = Pq + Spa 
and 9 ~ 9^ + S9 (pq is the stationary solution with 9'^ — 0) 
into Eqs. (JH) and keeping terms up to the first order in 
fluctuations {SpQ,S9), we have 



d_ 

d^ 

d_ 
d^ 



Spo=2XaPt,{l-pt,)59, 



S9 



-2 AXr, 



Uo d\a\ 



N dpo 



5P0: 



Po=Ph/ 



from which we find the small oscillation frequency lo as 
determined by the following equation 



If 



A + x|a|2)' 



c^2 = 4Aapg (1 - pg) 4Aa + 



Uod\a\ 



N dpo 



PO=PoJ 



In order to assure the dynamical stability of the system, 
up' should be positive. We find that in the region with 
three solutions, two of them are dynamically stable while 
the third one is dynamically unstable. This unstable 
state is shown by the the red dashed line in Fig. [2l it 
links the two stable ones, representing a typical example 
of bistability. 

In the region where the intracavity photon number is 
low, the interaction is dominated by the intrinsic s-wave 
scattering, which favors the ferromagnetic state in which 
Po = 0.5 for 771 = 0. In the region where the photon 
number is high, the cavity-induced effective Zeeman ef- 
fect takes a more prominent role which, for the choice 
of [/q < 0; favors a condensate in the rUg = ±1 mag- 
netic sublevels in which po becomes small. If a is fixed 
to a value independent of the atomic dynamics as in 
the case when it represents a strong off-resonant laser 
field [16], the system will experience a smooth crossover 
from the ferromagnetic interaction dominated phase to 
the Zeeman effect dominated phase as the strength of 
Uq is tuned. In our case, however, there is a first-order 
transition located within the bistable region as indicated 
in Fig. [2] This phase transition exists as a result of the 
cavity-mediated nonlinear atom-atom interaction. 



B. Non-equilibrium property: Critical Slowing 
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In this subsection, we study the spin-mixing dynam- 
ics of the system initially prepared in a state away from 
equilibrium. To begin with, we make use of Eq. ([S]) and 
rewrite Eq. Pa|) for to = as 



^) = 8A,po (1 - Po) [H-U (po)]-4 [H-U (po)]' , 

(6) 

where H is the energy of the system which is a constant 
determined by the initial condition. In the cavity-free 
model when U represents a constant quadratic Zeeman 
shift independent of po i Eq. ([6]) is known to support an- 
alytical solutions in the form of elliptic functions 12]. In 
our case, we have to resort to numerics to solve the above 
equation. As the system is conserved, the spin dynam- 
ics is expected to feature periodic population exchanges 
among different spin states, as in the cavity-free model 
with a homogeneous magnetic field [ij, 13 • 

Figure [3] shows how the oscillation period changes with 
cavity detuning (5c, where the period is obtained by solv- 
ing Eq. ([5]) numerically starting from the initial condi- 
tion (po = 0.1, 9 ~ O.IGtt) under the same set of param- 
eters that resulted in the equilibrium state in Fig. ([2]) 
with = 0. Here, cavity detuning 5c serves as a con- 
trol knob with which the departure between the initial 
non-equilibrium state (po = 0.1,61 = O.IGtt) and the clos- 
est equilibrium state (an in-phase state with 6 = Q) can 
be conveniently tuned. It plays a similar role as the mag- 
netic field in the study of spin dynamics in the presence 



FIG. 3: (Color online) Upper panel: Period of spin oscillations 
as a function of cavity-pump detuning 5c- Middle panel: The an- 
harmonic time evolution of po for the three peaks marked in the 
upper panel. Lower panel: From left to right, the phase-space con- 
tour plot of H corresponding, respectively, to the peak 1, 2 and 3 
marked in the upper panel. The black dots refer to the initial state 
of the system, while the white dots refer to dynamically unstable 
fixed points. 



of a homogeneous magnetic field. In the ferromagnetic 
case, it has been theoretically predicted [l2l that there is 
a single critical magnetic field around which oscillation 
period diverges. In contrast, the period as a function of 
5c in Fig. [3] exhibits three peaks around which the period 
(or the oscillation) exp eriences a dramatic enhancement 
(or slowing down) [39j. The spin population po as func- 
tions of time at three peaks are illustrated in Fig. |31 

To gain physical insights into these dynamics, we plot 
in the bottom of Fig. [3] the corresponding equal- iJ con- 
tour diagrams in the phase space defined by the conju- 
gate pair {9, Po). In a dissipationless system like ours, 
no matter how complicated the system dynamics may 
look in the time domain, it evolves along one such con- 
tour determined by the initial state (marked as a black 
dot in Fig. [3]). The critical slowing down takes place 
when the energy approaches a critical value He below 
which the contour changes its topology from an open to 
a closed line. In the pendulum analogy, it corresponds to 



the pendulum approaching the vertical upright position. 
The existence of a bistable region in our example makes 
the phenomenon of critical slowing down far richer. As 
can be seen, both the first and third peaks are located 
outside the bistable region, where only one attractor rep- 
resenting the stable state at 6' = exists, while the second 
one is inside the bistable region, where an unstable saddle 
point marked by a white dot coexists with two attractors 
a.t 9 — 0. Our results show that the oscillation period 
strongly depends on the cavity light field, the pump field 
can thus serve as a control knob for the spin-mixing dy- 
namics. 



IV. BEYOND SMA 

So far we have focused our discussion within the SMA. 
In this section, we will investigate the validity of the SMA 
and study the properties of the system when the SMA 
becomes invalid. 



Modulational Instability of a Homogeneous 
Condensate 
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In order to gain some physical insights into the valid- 
ity of the SMA, we first consider the case without the 
trapping potential and assume that the condensate in- 
side the cavity is homogeneous. In this case we have 
ipa = -y«Q exp (— i/iai — 1610), where the atomic density 
Ua now becomes position-independent, and the station- 
ary solution (n^, 6^) is still determined by Eqs. (|1]) 
at steady state except that Xa should be redefined as 

Xa = C2n/K. 

In order to check whether these homoge- 
neous states are stable against spatial modula- 
tion, we examine the Bogoliubov collective ex- 
citation spectrum by introducing small fluctua- 
tions around the steady-state solution. Inserting 
V'q = (y/n^ + S^^^a) e^p {-ifJ-it - i9^^) into Eqs. ([T]), 
where Sipa can be expanded in momentum space as 
(5-00 (r, t) = J2k ["a {t) exp {ik ■ r) + v^ {t) exp (-ik • r)] , 
we obtain a matrix equation idx/dt = A^x for vector 

X ~ (u+,Mo,'u_, w^, Wq, w!l) where 7W is a matrix 
given in the Appendix. The Bogoliubov modes are 
then given by the eigenvalue equations A^x = fiwx, 
where w represents the excitation frequency if it is real 
and signals modulational instability with a growth rate 
Im (oj) if it is complex. 

In the absence of cavity, the homogeneous ground state 
of a ferromagnetic ^^Rb condensate is stable against spa - 
tial modulation even with a finite magnetic field [20|. 
This, as we shall show, will not be the case when the 
cavity is introduced. To illustrate this, we numerically 
diagonalize M to investigate the properties of Bogoli- 
ubov excitations. For the spin-1 system as we considered 
here, there will be three branches of Bogoliubov excita- 
tions — two gapless branches and one gapped branch. 



FIG. 4: (color online) (a) The grey area shows the range of 
wavevector k corresponding to the unstable excitations of the self- 
consistent ground-state of a homogeneous rubidium condensate- 
cavity interacting system. The red lines refer to those with the 
maximum instablity growth rate, (b) Spin domain width versus 
cavity-pump detuning Sc- 



Our numerical calculations reveal that one of the gapless 
branches will become unstable for certain values of the 
wavevector k. Figure Sl^a) shows the range of unstable 
excitations associated with \k\ £ [0,fcin], and those with 
the maximum instability growth rate are represented by 
the red lines. Furthermore, the Bogoliubov eigenvectors 
of these most unstable modes are found to take the fol- 
lowing form 



,wj ex (-0.5,0,0.5) or (0.5,0,-0.5), 



which describe the spin waves with spin angular momen- 
tum ±h. The exponential growth of these modes tends to 
induce spontaneous magnetization, and spin domain will 
be formed as a result of the competition between local 
spontaneous magnetization and the conservation of the 
total magnetization. The size of the spin domain may 
be estimated by the inverse of the wavenumber 2n/k,^, 
which is plotted in Fig.|4|b). 

It is important to note that if the total size of the con- 
densate is small compared to the domain width estimated 
above, the instability will be suppressed. 



B. Spin Domain Structure 

Equipped with the insights gained from the study 
of a homogeneous condensate in the previous subsec- 
tion, we are now in the position to explore the effect 
of cavity-induced atom-atom interaction on spin-domain 
formation in a trapped condensate. For simphcity, we 
consider a cigar-shaped trap with a harmonic trap po- 
tential Vt (r) = nia [uj"^ {x^ + J/^) + ^l^^] /2 in which 
the transverse trap frequency ujj_ is much higher than 
the longitudinal trap frequency uJz- This allows us to 
introduce a longitudinal wavefunction (l)a{z,t) via the 
ansatz ipa (r, t) = (J)a_ {x, y) 4>a [z, t) exp {—2iu]±t), assum- 
ing that the transverse wavefunction 0_l (cc, y) always 
remains in the ground state of the transverse poten- 
tial. Following the standard approach (see, for example, 
Ref. [Ull^), we simplify Eqs. dTa]), ^ and © into a 
set of equations for (pa {z, t) 



ih(j)±^ C + Uo\a\ + C2 {p± + po - Pt) 

ih(po^ C + C2{p++ P-) 0o + 2c20+(/)- 



(7a) 
(7b) 



which describe an effective ID trapped system, where 



£ 



2ma dz 



2 ' -l^^z^ +cqp, 



with Pa = \(t>a\ , p = p+ + po + P-, and Co(2) 

Co(2)TOaW±/27rft. 




FIG. 5: The ground denstiy profile of a Rb condensate trapped 
in a unidirectional ring cavity. Here the distance z is scaled with 
Qs = y'h/maiOz, and the parameters used are specified in the 
main text. 



In our calculation, we set the trap frequencies as 
uj^ = (27r) 240 Hz and uj^ = (27r) 24 Hz [ll], and other 
parameters same as before. The Thomas-Fermi radius in 
the z-direction is then about 24 /im,. In the numerical 
simulation, we obtain the ground state in a self-consistent 
manner by propagating Eqs. ([7]) in imaginary time sub- 
ject to the constraints set by the conservation of both the 
total particle number and the magnetization. The results 
are shown in Fig. [5j From the numerical simulation we 
find there exists a critical value of the cavity-pump de- 
tuning Ac ~ —3.5. The ground state exhibits a typical 
spin domain structure when Sc > Ac, in which mg = 1 
and —1 states occupy the opposite ends of the longtitu- 
dial trap. While for Sc < Ac, all three spin components 
are completely miscible with no spin domains forming, 
and the ground state is well described by the SMA. This 
is a clear proof that the cavity light field can be used to 
control spin domain formation in the condensate. The 
mechanism lies in the fact that the domain width can be 
significantly modified by tuning the cavity-pump detun- 
ing ^c, as we have shown in Fig. |H[b). When 6c > Ac, 
the domain width (around 14 fj,m) is smaller than the 
condensate size and spin domain can be formed. While 
for Sc < Ac, the domain width is larger than the size of 
the condensate, then the domain formation instability is 
suppressed. 

At this point, we comment that spin domains was first 
observed in the ground-state of a ^'^Na antiferromagnetic 
condensate in the presence of the magnetic field gradient 
Pi . Later studies }17nl9i] discovered that a ferromagnetic 
spinor condensate initially prepared in an excited state 
will be subject to dynamical instability and lead to spin 
domain formation, while antiferromagnetic ones are dy- 
namically stable. The experiment of Ref. 2] displayed 
the spin domains formation in a quenched ^^Rb ferro- 
magnetic condensate. Recent work [20[ clarified that for 
a spin-1 condensate subject to a homogeneous magnetic 
field, the ground state exhibits domain formation only 
in antiferromagnetic condensates, but not in the ferro- 
magnetic ones. The significance of our work here is that 
spin domain structures can also be created in the ground 
state of a ferromagnetic condensate with the aid of a cav- 
ity. This can be traced to the effective spin-dependent 
atom-atom interation induced by the cavity. 



V. CONCLUSION 

In conclusion, we have studied the mutual interaction 
of a ferromagnetic spin-1 condensate with a single-mode 
cavity. The intracavity light field and condensate wave- 
functions are calculated self-consistently. The cavity- 
mediated effective interaction gives rise to a variety of 
new spin-depedent phenomena. Under the SMA, both 
the equilibrium properties and non-equilibrium dynamics 
are investigated in detail. We show that the system can 
display bistable behavior. By tuning the cavity-pump 
detuning, the spin-mixing dynamics can be manipulated. 



We also discussed the situation when the SMA becomes 
invaUd, and found that phase transition among different 
spin components can occur in the ground state which 
leads to spin domain structure. All these effects can be 
readily tested in experimens. The cavity-spinor conden- 
sate system can provide a new platform for the study 
of cavity nonlinear optics and the properties of spinor 
condensates. 
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Inserting Tpa = {^/n^ + 5ij:a)eyi]i{—i^%t — i9%) and 
a = a* + (5a into Eqs. ([Ta| and (jlbl) where a* is 
the steady-state value of Eq. ^ corresponding to the 
equilibrium solutions {n^,0^), in the homogeneous case 
{Vt = 0), keeping terms up to first order in 6ipa and 5a, 
we obtain 



J 



iMV'i = -h'^\7'^/2ma - /i± + Uo |a1^ + 2 (cq + C2)n% + (co + cj) ng + (cq - ca) n^ St/j.^ 

+ [{cq + C2) nl + C2ng exp {-ie")] 5^^ + (cq + ca) i/^g'^i ('^^o + ^V'S) 
+ (co - C2) y^?l4 {5i>^ + 5i>*^) + 2c2 v/^Ig^f^V'o exp {~i6') 
+ Uo^/n^[a''5a* +a''*5a), 



(Al) 



zMV'o = [-/i^V^/2ma - Mo + (^o + C2) {n\_ + nl) + ^cmt)] 5ipo + con'„5Tpo 
+ (co + C2) [^^4^ {5^J+ + 5i>+) + V^^ {5ip^ + 5iP*_)] 
+ 2c2 ( v/^4^(5Vo + V'^^ng'^V^- + VnS^5tl;+) exp {iO") , 
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(A2) 
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(A3) 



where V = N/n is the volumn of the condensate and the 
use of Eq. ([2]) has been made in arriving at Eq. (|A3p . 
Finally, by combining Eqs. (KJl . (|X2|) and (|X3)) . we can 
construct matrix A^ in a straightforward way. 
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